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Renyi and von Neumann entropies quantifying the amount of entanglement in ground states of 
critical spin chains are known to satisfy a universal law which is given by the Conformal Field 
Theory (CFT) describing their scaling regime. This law can be generalized to excitations described 
by primary fields in CFT, as was done in reference []]], of which this work is a completion. An 
alternative derivation is presented, together with numerical verifications of our results in different 
models belonging to the c — 1,1/2 universality classes. Oscillations of the Renyi entropy in excited 
states and descendant fields are also discussed. 

I. INTRODUCTION 

Entanglement is a fundamental property of quantum states according to which measurements in one of two sub- 
parts of a quantum system can (immediately) condition the state of the other part, which may be far away in space. 
This property was at the origin of Einstein's criticisms to quantum mechanics in his EPR article in 1935. Today, 
entanglement is one of the most studied topics in physics and many research lines are related to this concept. It is the 
property that makes possible quantum computation, teleportation and quantum information processing [2J. On the 
other hand, the irruption of quantum information methods and concepts in the condensed matter physics community 
has led to one of the most fruitful research areas of the last decade [|| . Measurements of entanglement between spatial 
or more abstract parts of many-body systems and statistical models in their ground state have been proved to unveil 
essential information characterizing their phases. After all, quantum correlations, or entanglement, in ground states 
(gs) of many-body systems are responsible for the onset of coherent phases of matter at zero temperature such as 
superconductivity and Hall states. Probably the most studied quantity in this context is the amount of entanglement 
between two spatially separated parts of an extended system. Calling A one of these parts, the reduced density 
matrix, pa, of the ground state under study, is obtained by tracing the degrees of freedom of the complementary of 
A. If A results in a mixed (pure) state, then A and its complementary are said to be in an entangled (product) state, 
the amount of entanglement being normally quantified via the von Neumann entropy of pa- 

Si{A) = -tr p A lnpA, (1) 

called entanglement entropy. Alternatively, the n-th order Renyi entropy is used, n > 1: S n (A) — :rf— lntr^ pA n , 
which contains information on the spectrum of pa, the entanglement entropy being lim„_j.i S n . The interest of these 
quantities is threefold. Firstly because of the area law satisfied by ^ (see references in [J]), stated in its first version 
by Bombelli et al in 1986 and by Srednicki in 1993 for a free scalar field, and related to the problem of information 
loss in black holes 5] , [6] . Roughly speaking, it states that ground states of Hamiltonians with short-range interactions 
are such that their entanglement entropy is proportional to the area of the hyper-surface separating A from the rest 
of the system, hence proportional to £ in d dimensions, if I is the linear size of A. This important property serves 
to characterize the very small fraction of the Hilbert space accessibly to "physical" ground states. In one spatial 
dimension (ID), the area law implies that the entropy of massive ground states is bounded and independent of the 
length I of the linear part, or block, A. Being the entanglement entropy proportional to the minimum amount of 
information needed to describe the partition A, this was used (see references in [3]) to interpret the efficiency of the 
Density Matrix Renormalization Group (DMRG) algorithm in massive phases Q, and the lack of this efficiency when 
dealing with gapless phases, a case in which the area law is to be corrected with a logarithmic term in one dimension, 
leading to an unbounded entropy in infinite systems. In the same spirit, measurements of entanglement are essential 
in the development of novel [7[ efficient algorithms for the simulation of quantum many-body systems in one and 
two dimensions. Entanglement measurements, from a broader perspective, may provide relevant information about 
the physics of the system under study. Just to give an example that will be cited in the body of the article, the 
oscillations of the Renyi entropy in critical ID chains have been proved [9j to be related to the Luttinger parameter 
entering in their continuum description, which, in principle, could be inferred in this way through quantum information 
measurements in particular finite-size models. 

Finally, the interest of entanglement entropy comes also from the fact that the mentioned logarithmic behaviour 
of the entanglement in ID has been proved to be a universal property of critical systems, captured, among other 
properties of entanglement, by their underlying Conformal Field Theory (CFT). Since the work of Polyakov et. 
al. [Io| . conformal symmetry of critical two-dimensional systems has been exploited to infer the universal form 
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of correlators, finite size scaling of energy and momentum, critical exponents and several properties of stochastic 
evolution of interfaces [n], EH- Moreover, in 1994, Holzhey et. at, and Calabrese et. al. in 2004 showed in a CFT 
context that, in a system of length N and periodic boundary conditions (PBCs), the Renyi entropy takes the universal 
form 0,(11], [3: 
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where c is the central charge of the CFT describing the system scaling limit and 7„ is a non-universal constant. This 
statement was generalized to account not only for the ground state entropy but also for low-energy excitations in 
reference [l[ , a generalization which supposes a further prediction of CFT in the held of quantum critical phenomena. 
By conformal invariance, the hnite-size spectrum at criticality has a universal structure determined by the Virasoro 
algebra satisfied by the generators of the conformal transformations. Each state in the low-energy spectrum exhibits an 
energy and momentum finite-size scaling determined by the conformal weights, or, roughly speaking, the eigenvalues 
of the Virasoro operators of the given state. In the same spirit, in reference [l[ it was proved that the Renyi entropy 
of excitations defined by primary fields are universally related to conformal properties of the operator defining the 
targeted excitation. In particular, consider the quantity: 



4 n) =exp[(l-rO(^(x)-Sf(z))], (3) 

defined in reference It quantifies the excess of entanglement of a state |T) related to the primary field T, with 
respect to the gs. This quantity [1] related to the 2n-point correlator of the operator T and its conjugate in the 
cylinder: 



(T(0)Tt(27ra))S y 



F<t n \ x ) = n- 2n( > h +v Ulj = /; 2 v » //cy . (4) 



This is the main result we will alternatively derive in the following section, and that will be numerically checked in 
subsequent sections for different models. 

Entanglement of excited states has been also considered in [l6[ , where universal scaling of the negativity of excita- 
tions in the XX Z critical chain was shown. In [l7j it was shown that a violation of the area law should be expected 
for the low lying excited states of critical chains, and in [l8[ it was considered the entanglement of very large energy 
excitations in the XY and XX Z models. In [19], [2(|, the law (j4j was applied to study systems with continuous 
degrees of freedom, where it turned to be an accurate description of the finite-size system entanglement. 

The aim of this work is to complete the results of reference [1]. In section [TT| an overview of entanglement of 
excitations in the XX model is given. Section [IIII is a derivation of our results in the framework of the Calabrese and 



Cardy computation [15[ (in [1|, the approach of [13( was used instead). Section [TV] presents some of our results for the 
bosonic theory and some numerical illustrations for three models in this universality class: the XX (free fermions), 
XXZ and excluded-volume-AA spin chains. The section concludes with a study of the Reny entropy oscillations in 
this theory. Section [V] is a study of the free Majorana CFT with numerical tests in the quantum critical Ising chain. 
A discussion about the entropy of descendant fields is given in section PVTl 



II. A FIRST SIGHT ON ENTANGLEMENT OF EXCITATIONS IN THE XX MODEL 

As a first illustration of the behaviour exhibited by block entanglement in excited states of spin chains, let us 
consider the XX model for N interacting spins occupying the positions of an TV-site lattice and whose interaction is 
defined by the Hamiltonian: 



"w ^7»;<i • «.,)• (5) 

j'=i 

where o~j are the spin Pauli matrices acting on the j-th spin. This well-known model [2l| is a paradigm of spin chain 
describing quantum magnetism, it also emerges as the strong on-site repulsion limit of the boson Hubbard model. 
As discussed in appendix [S] and in reference [2l[ , the XX model is integrable and it can be mapped, through a 
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Jordan- Wigner transformation, into a problem of lattice free fermions. According to such a mapping, the spectrum 
of the Hamiltonian ([5]) coincides with the spectrum of a free-fermionic Hamiltonian: 



H = Ej djdj + constant (6) 

j' en 

where {dj,dj,} — 8j t j>, are the creation and destruction operators of fermions with momentum kj — 2%j/N, Ej = 
— cos kj is the free-fermionic dispersion relation, and fHs a set of N integers or half- integers such that the resulting 
momenta kj are the ones determined by the periodicity or anti-periodicity of the boundary conditions in the fermionic 
problem. In other words: each one of the 2 N eigenstates of is associated with one of the 2 N fermionic eigenstates 
of©: 



n F <N 



n 



(7) 



(being |0) the fermionic vacuum annihilated by the d's and {m,,}™^ = K C il, one of the 2 N subsets of O), and the 
association is such that the eigenvector in correspondence with the momentum set K has energy eigenvalue: E — 
— X^gk cos explained in the appendix, the Jordan- Wigner transformation is such that a fermion Hamiltonian 

with anti-periodic boundary conditions (APBCs) is related to the Hamiltonian (J5J) with periodic boundary conditions 
(PBCs) for even uf, and with APBCs for odd uf- We will impose APBCs to the fermions. The corresponding O is: 



In this framework, one has that the ground state of ([3]) corresponds to the Fermi state of © with hf = N/2 
fermions: 



\n F )= J] 4^.|0). (9) 

0<j < (« F -l)/2 

From now on we will denote \m) the Fermi state of m fermions, as in equation ([9]) (even m is supposed). Let us 
now consider some low-energy excitations in this model. Due to the conformal invariance, the low-energy excitations 
present an excess of energy with respect to the ground state which is: 2nA/N, being A the conformal dimension 
of the excitation. The excess of energy vanishes in the large- N limit: the model is critical. As a critical model, its 
ground state \N/2) satisfies the law ©, one can ask whether the entropy of low-energy excitations coincides with the 
gs entropy in the thermodynamic limit, as happens for the energy. It turns out that there is a class of low-energy 
excitations for which this is indeed true. This type of excitations will be called compact since, as the ground state, 
they do not exhibit holes in momentum space or, in other words, for them the set K is composed of consecutive 
momentum quanta differing by 2n/N. Otherwise the excitations will be called non-compact, and they are such that 
their entropy is larger than that of the ground state. In what follows we illustrate with examples how excitations in 
the compact class presents an entropy equal to the gs entropy, up to the oscillations described in Q, with the same 
non- universal constant 7„. Afterwards, some non-compact states will be analyzed. These results will be justified with 
CFT arguments in section HVl 

As a first example of compact states we consider the excitation obtained by removing the highest momentum fermion 
in the XX model, or d nF /2-i/2\ n F >, with np = N/2. In figure [T] we present the n = 2, 3,4-Renyi entropy of such 
a state (which is labelled as (a)), together with the ground state entropy. The different behaviour of the oscillations 
(which are present only for n > 1, and absent in the N oo limit) is the only difference between ground and excited 
state entropies. As a further illustration, we present the state obtained adding a fermion below the left Fermi point 
and another one above the right Fermi point (i.e., the state ^ F / 2 +i/2 ^-n F /2-i/2l ni? ))- This state is labelled as (b) in 
figured] A final example of compact states is provided by the Umklapp excitation (c.f. table [J), consisting in moving 
the fermion with most negative momentum to the right of the positive Fermi point, i.e. to eL(„ F _i)/2 < ^(n J r+i)/2 1 77 -^) ■ 
This can be proved[l| to have exactly the same entropy S n as the ground state for all values of n (with exactly the 
same excitations), since such a shift in momentum space k — > k + 2n/N amounts to a phase shift of the wave-function 
in position space, leaving the reduced density matrix unchanged. 
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name of excitation 


field 


(h,h) 


state (n F = N/2) 


(h:p) 


N — 8 example 


ground state 

(a) 

(b) 

Umklapp 


1 

e i<t>-i<t> 


(0,0) 
(1/2,0) 
(1/2,1/2) 
(1/2,1/2) 


\n F ) 

d(n F -l)/2\nF) 
^n F /2+l/2 ^-np/2-1/2 ™ F ) 

d-(n F -i)/2d\ nF+1)/2 \n F ) 


(0 

(10 


o o • • • • oo 

o o • • • o oo 

o o o • • • «o 


particle-hole 
R-L particle-hole 


idcj) 


(1,0) 
(1,1) 


d (n F -i)/2 d ( +iy2 \n F ) 

d ( n p-l)/2 d („ F + l )/2 
<*-(n F -l)/2<*l ( „ F + i )/2 l™F> 

d(n F -i)/2d'\ nF+:i)/2 \n F ) 


(1:1) 
(1:2) 


o o a a a o ao 
o a o a a o ao 
o o a a a o oa 



TABLE I: A summary of the mentioned excitations. The horizontal line separates the compact states from the non-compact ones. The 
notation (h,p) applies only for chiral excitations, and the corresponding conformal fields are shown for primary states only. 

For the sake of clarity let us introduce the following notation: [hi hi ■ ■■ ■ P1P2 ■■■) refers to a chiral excitation 
with holes in the hj's allowed momentum values below the right Fermi point and particles in the pj's momentum 
values above it, i.e., to the state: 

nd \ F -± , dnp+i, \n F ), 
— % hPj 2 n * 

j,k 

in such a way that the ground state is denoted as ( : ) and the (a) excitation is (1 : ). 
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FIG. 1: Renyi entropy S r , r = 2,3,4 for three states in the XX model with N = 200 sites. The states (a) and (b) are described in the 
text and in table U The continuous curves are the universal function J2}- 

We will now focus in non-compact excitations, such that holes are created in momentum space. These states do 
not obey the law ([2]). The lowest excitation of the XX model with non-zero momentum has excess of energy and 
momentum equal to 2n/N, and in the fermionic language it corresponds to the destruction of the fermion below the 
Fermi point and the creation of a fermion immediately above it. In other words, it is the (1 : 1) state. We will call 
this excitation a particle-hole excitation. This state is found to exhibit excess of entanglement AS(£) — S(£) — S gs (£) 
which is larger than zero. Figure [2] represents the l/N « 1 regime of the excess of entanglement entropy for such 
a state, together with other particle-hole like excitations. Since the works [l3l - [l5l ]. one knows that, for small l/N, 
the ground state of a critical system satisfies S\(t) ~ | \nl. For excited states in the XX model, one observes a 
correction of the type [l|: 

2 _2 / p \ 2 

ASi (I) = — ( h + ~ h )[x) + ° [(£/A ° 4] (10) 



ground state □ 
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where ASi is the excess of entropy with respect to the ground state and h, h are the conformal weights of the operator 
corresponding to the excitation (equation (|10[) will be justified in section Hv]) . As a further example consider the state 
d(n F -i)/2 ^(n F -i)/2+ml ni? )' sucn that the fermion nearest to the Fermi point has been displaced m momentum quanta 
to the right. On inset of figure El one observes for this state how curves for different m collapse when the quantity 
ASi/m is plotted, confirming also in this case the law (TTU)) with h = m, h = 0. 

It is worth to stress that although some of the states studied in figure [2] present the same low- a; dependence of their 
entanglement and Renyi entropies, their entropy do not coincide for larger values of x in general. This is illustrated 
in figure El in which we show the function defined in ([3]) for the excitations studied in figure [2] 

In section HVl we will show how some of the results presented in this section can be justified with CFT arguments. 
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III. CFT APPROACH TO THE ENTROPY OF PRIMARY FIELDS 

In reference [l[ we obtained a general formula for the Renyi entropy of excited states in CFT following the approach 
of Holzhey, Larsen and Wilczek (HLW) who computed the entanglement entropy of the ground state of CFT systems 
with periodic boundary conditions. The HLW approach was later generalized by Cardy and Calabrese (CC) to deal 
with more general situations, as finite temperature, open boundary conditions and several disjoint intervals. It is thus 
of great interest to re-derive the result of reference [l[ using the CC approach. This is the aim of this section. 

Let us start by briefly recalling the main steps of the CC approach. The general set up of the problem is a lattice 
1+1 quantum field theory with local commuting observables {<f>{x)} with eigenvalues {4>(x)} and Hamiltonian H. The 
thermal state p at inverse temperature (3 has matrix elements 

({<P"(x")}\p\{<f>'(x')}) = ^({0'V')}|e-^| W)}) (11) 

where 

Z(J5) = Tre~ pH (12) 
is the partition function. Equation (fTTj) can be expressed in terms of an euclidean path integral as 



{{<t>"{x")}\ P \{4>'{x')}) =y^\ ^ r )i n wo*, -pm - n pi 2 ) - <w)) e ~ SE 

where Sb = S-p/i ^ T with L e the euclidean Lagrangian. The normalization factor in ([lip guarantees that Tr p = 1 . 
The partition function Z(j3) is obtained doing the path integral with the identification <f>'{x) = <j>"(x) at r = — /3/2 
and r = (3/2, and integrating over these variables. The geometry of the integration surface is a cylinder of length (3. 

Let us now take a system A given by the interval (u, v). One wants to compute the reduced density matrix p{A) 
by tracing over the points not in A, that we call B. This operation amounts to gluing those points in (|13|) and doing 
the integral over them. The effect is to leave a cut (u,v) along the line r = —(3/2. To compute tn p n (A) one uses the 
replica trick. One first makes n copies of (fL3|) . labelled by k = 0, . . . , n — 1 and glue them together cyclically 

<t>' k {x) = <j>'l +l {x), (fc = 0,...,n-2), c(>' n _ x {x) = <%{x), VxeA (14) 
The path integral on this n-sheeted geometry is denoted as Z n (A) and hence 



trp «(^) = M4) (15) 

CC argue that the LHS of ([i"5j) is analytic for all 5Rn > 1 and that its derivative respect to n in the limit n —> 1 + 
gives the entropy 



Sa = — lim — — "1. - (16) 

n-»-i dn Z((3) n x ' 



Alternatively, one can define the Renyi entropy S n (A) for the interval A, 



S n (A) = -^—logtrp n (A) (17) 



and compute the entanglement entropy as 



Sa = lim S n (A). 

n— >1 



(18) 
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Let us now turn to the excited states in CFT that correspond to primary states. If |0) is the vacuum state then 
the incoming state generated by a primary operator T(z, z), with conformal weights h and h, is given by 



|T)= Urn T(z,z)|0) (19) 

z,z— >0 

while the outgoing state is given by 



(XI = lim z- 2h z- 2h (0\r^ [-,-) (20) 

z,z^0 \z Z J 

where is the operator conjugated toT, i.e., TxT^=I-| — •. In the latter two equations we have used the radial 
quantization so that z = and z — oo correspond to the infinite past and infinite future respectively. The primary 
operator corresponding to the vacuum |0) is, of course, the identity, i.e. T = I. The density matrix of the ground 
state is given by the limit ft — > oo of equation ([13")) , which we write as 



({<P"(x")}\pi\W(x')}) = - / [d<f>(x, r)] [] S(4,(x, -oo) - 4>'{x')) J] 6(<f>(x, oo) - <j>"{x")) e' 3 - (21) 

•' X X 

where Z — lim^oo Z(fi) and the action Se is computed over a cylinder of infinite length. Based on equations (|19[) 
and (|20p we find the following expression for the density matrix of the primary state |T), 



({</>" (x")}\pr\W(x')}) - C / [d<j>(x, r)] J] 5(ct>(x, -oo) - <t>'{x')) J[ 8(<j>(x, oo) - cj>"(x"))T(0, -oo) T*(0, oo) e~ s - 

J X X 

(22) 

which reduces to (|2ip for T = I. In this equation we have parametrized the primary fields in terms of the space-time 
coordinates, i.e. T(tr, r), etc. not with the radial coordinates, i.e. T(z,z), as in equations (IT9l) and (|20|) . but there is 
a one-to-one correspondences between the two. C is a constant which will be fixed below. Following the same steps 
as for the derivation of equation (fT5|) . one finds that the trace of the reduced density matrix p^{A) is given by 



tr ft{A) = C"Z„(A)(T (0, -oo) Tj(0, oo) . . . T n _i(0, -oo) T^O, oo)) K „ (23) 

where lZ n is the n-sheeted Riemann surface that results from the sewing of the n copies of the original cylinder along 
the cuts associated to the interval A. The label k = 0, ... ,n — 1, of the primary field Y^, and its conjugate Tj,, 
denotes the sheet they belong to. The constant C can now be fixed imposing the normalization of the reduced density 
matrix, 

Trpr(A) = C r Z(T (0,-oo)Tj(0,oo)) TCl = l (24) 
Eliminating C from this equation and plugging its value into (|23p yields 



+ Z n (A) (T o (0, -oo) 4(0, oo) . . . -oo) 1^(0, oo))^ 

^'"Z 5 (T (0,-oo)T 5 (0,oo))^ (25) 



Finally, using equation (1151) in the limit (3 — > oo, one gets 



F („) _ trp^(A) = (T o (0, -oo) Tl(0, oo) . . . T n _!(0, -oo) ,,(0, oo))^ 
^Pi(A) (T (0,-cx>)TS(0,oo))^ 

Hence, the ratio of the trace of the reduced density matrix of a primary state and that of the ground state, is given 
essentially by a 2n point correlator of the primary field T and its conjugate field on a n— sheeted Riemann surface. 
In the previous derivation we have not made used of the fact that A is given by a single interval (u,v), so equation 
(1261) also holds when A consists in the union of multiple intervals. 
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FIG. 4: The cylinder w with the cut A is mapped via l|27|l to the center panel (we took x = \u — v\/L = 1/8). The transformation l|30|l 
with n = 3 has the effect shown in the right panel. The red crosses indicate the images £oo, f^, Zfc^, z' k 3 of the points ?«oo, w'^. 



Let us next compute the correlators appearing in equation (|26p in the case where A = (u,v). To this end we shall 
first introduce the complex coordinate w = a + ir, which parametrizes an infinite cylinder of length L, i.e. a = a + L. 
The interval A will be identified with the domain u < w < v. This cylinder can be mapped into the complex plane 
by the conformal transformation 



c = 



sin 



7r(w — v) 



(27) 



The points w — u and w — v are mapped into £ = and £ = — oo respectively, so that the interval A corresponds 
to the negative real axis < 0. Moreover, the regions 3w ^ are mapped into the regions ^ 0. The points 
(0, — oo) and (0, oo), appearing in equation (|26[) correspond to the coordinates Woo — —zoo and w'^ = zoo, and, by 
equation l|27p. to the ^-coordinates (see fig HI) 



Coo ^ 



So 



t = v- 



(28) 



Notice that x belong to the interval (0, 1) since I € (0, L). The Riemann surface 1Z n is constructed out by gluing 
cyclically n copies of the complex £-plane along the cut A. The uniforming parameter associated to lZ n is given by 



l/n 



' . n(w-u) \ 1 / r ' 

sin — ^-r — - \ 



(29) 



Consequently, each of the points Coo and C£o give rises to n points in the Riemann surface lZ n with coordinates 



~{x+2k) 



k = 0, . . . , n — 1 



To compute (|26p we first write it as 

f!^\a) = f^\x) 



lim 

■w— > — ioc 



(Uk=o T k(w, w) T\(-w, -w)) 



(T (w,w)tJ(-w,-w;))™ i 



(30) 



(31) 



The 2n point correlator on lZ n , appearing in the numerator, can be transformed into a correlator on the complex 
z-plane by means on the conformal transformation of the primary field T 



T(w, w) 

and a similar equation for TK From equation (|29p one finds 



—\ ( —\ X(z z) 
dw J \ dw J 



(32) 
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dw n L 
which in the limit w — > ±ioo becomes 



"' ; 4lr sm(nx) ( e ¥(--«) _ e -¥(--«))-i( e ¥(— ) „ e -¥(---))-i (33) 



where 



dz Zfc n A / , \ 

w — > —too, z — >• Zfc „ — > A (34) 

aw n 

I dz z k.n t /or N 

w —> too, z — > z kn => — > — — A (35) 

dw n 



A-rr 

A = — am(nx)e' 27rlwl/L e t7r( - u+v ' >/L (36) 



In this limit the conformal transformation (|32[) yields 



W -t -tOO, Z -> Zk, n T k (w, w) -> ( A) A) T(zk, n ,Zk,n) (37) 

h / —1 \ h 



w^ioo, z 2 fen =J> T£ («>,«/) ->• I A J AJ T T (z fe ,„,2; fc! „) (38) 

Plugging these equations into ([31]) , the factors proportional to A and A cancel out and one is left with 

p(n)i \ _ -2n(h+h) IIfc=0 z k,n^\,n Z k,n^k,n (TTfc=0 ^ ( z k,n, z k.n) ^H Z k,n> Z k,n))c , . 



/'-ft g ft ,/ii t ^\" (T(z ,i'^o,i) Tt (4,i)^o,i))c 

^ Z 0,1 Z 0,1 0,1 0,1 J ' ' ^ 



To simplify this expression, we can make a further conformal transformation from the complex plane to a cylinder 
of length 2tt, 

z = e lt -> T(t, £) = ^i h ~~ h ) z h z h T(z, z) (40) 

which allow us to write (15"9")) as 



p(n)/ \ _ _2n(ft+ft) (Ilfc^Q T(tk,n,tk,n) T^f^, 4,rJ) 

T (T(to 1 i,<o,i)Tt(* / o, lJ *o,i)) n 



cy 



(41) 



where 



t k , n = -(x + 2k), t' kn = -(-x + 2k), fc = 0,l,...,n-l (42) 
Since the latter parameters are all real we may as well write (|39p in a compact form as 

M M - Mh+ H) (n'£oT(^(.x + 2fc))Tt(iL(-x + 2fc))) cy 

where T(t) = T(t,t). After a shifting of the points in which the correlator is evaluated, this equation become (j4j, 
which was derived in reference [l| using the HLW approach. The function F r (x) satisfies the following identities 



" } (x) = (f^ (xj) * = 4f) (*) = (1 - x) (44) 

that can be derived from equation (|39|) using the transformation properties of correlators under inversion z — > —1/z 
and rotations z — > e %a z. The last equality in (|44p reflects the well know fact that the Renyi entropy S n (A) is the same 
as S n (B), where B is the complement of A, which amounts to the replacement I — > L — i, that is x — > 1 — x. 
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IV. BOSONIC THEORY: XX, XXZ AND EXCLUDED-XXZ MODELS 

The first theory under study is a massless bosonic field ip(z, z), with action: 



•A-W\ = J dz dz d z tp dzip 



(45) 



This is a CFT with central charge c = 1 and two types of primary fields. The first type is given by the vertex 
operators 



^i(a<t>+acf>) . 



(46) 



where <j> , <f> are the chiral and anti-chiral parts of the bosonic field, solution of the equation of motion: ip(z, z) = 
4>(z) + 4>(z). In the last equation : • : denotes normal ordering and a, a are real numbers. The vertex operator has 
conformal weights (h, h)=(a 2 /2, a 2 /2) 



(Va 1 ,a 1 (*l,«l)Va aj a a (2&,22)> = ( z i ~Z2) aia2 (21 -z 2 f l£ 



(47) 



where (•) denotes the correlator in the complex plane, and the neutrality condition implies that the correlator is 
different from zero only if oti + cti =0 and the same for the a's. This result is indeed more general. Considering 
holomorphic fields a = only, one has {zij = Zi — zj): 



(V ai (z 1 )...V an (z n )) = l[[z ij r^ 



(48) 



i<3 



if otj = 0, and zero otherwise. In a cylinder of length 2ir parametrized by w = — i In 2 this correlator is: 



(V ai (w 1 )...V a Jw n ))cy = II[ 2 

sin(wij/2)] c 



(49) 



l<3 



evaluating ([49]) in the points indicated in (|4]), i.e., in the set {2irj/n, 2ir(j + x)/n}^ , j = 0, . . . , n — 1, one has : 

sin(7rm/n) 2 



This function is constantly equal to one: 





sin(7ra;) 


n n—l 




n 

rn—1 


sin(7ra;/n) 



7t(m— x) • ir(m-\-x) 
sin — 1 L Sill — i 



(50) 



F^(x) = l Vn,a, 



(51) 



a fact that can be proved by comparing the zeroes in x of the numerator and in the denominator of (1501) . There are 
zeroes of order n in x G Z in both the numerator and the denominator. As analytic functions, this implies they are 
proportional, the proportionality constant being 1 as can be immediately seen in the x — > limit. We thus observe 
that the entropy of excitations corresponding to vertex operators coincide with the ground state entropy according to 
the CFT prediction (gj). 

Let us now focus our attention in the other primary field in the theory: id<f>, with conformal weights (h, h) = (1, 0). 
We will use the correlators of idtfi in the plane: 



1 



(d<l>(zx)d<t>{Z2)) = —z 



212' 



(52) 



and the Wick theorem: 



2n 



(11 idftzi)) = Hf 

3=0 



(53) 



l<i,j,<2n 
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n F = 20 N = 120 
n F = 250 N = 500 
riF = 250 N = 500 
eft predictions 




x = £/N 

FIG. 5: Renyi entropy ratio of the particle-hole excitation in the XX /free fermion model, n = 2,3, and two different sizes, vs. the CFT 
prediction equations H571 1581) (continuous lines). For n = 2, a system with filling fraction 1/6 is also shown. The inset is a zoom of the 
region defined by the small rectangle over the data-set in the body of the figure. 



where Hf [] denotes the Haffnian (S^n is the permutation group of the 2n indexes): 



1 

Hf ( A ) = II A 7r(2i-l),7r(2i), 



(54) 



and where the matrix in ()53[) has a null diagonal. This matrix being dependent only on i — j, can actually be written 
as a determinant 1231: 



Hf 



1 




' 1 " 




= dot 


.**3 _ 







(55) 



In the cylinder w — — zlnz, and taking again the 2n coordinates Wj as in (|4]) one gets: 



F$(x) = n~ 2n [sinTrx] 2 ™ det 
The n — 2,3 values of are reported: 



sin (wij/2) 



l<i,j<2n 



(56) 



F^(x) = l-2v z + 3r 4 - 2r b + r 



6 , _8 



v 7 729 



l-4s 



2s 3 



2s 3 



3s 2 



(57) 
(58) 



where r = sin7nr/2, s = sin7rx/3, s± = sin7r(l ±x)/3. These values will be compared below with the numerical 
values of obtained for the finite-size states \id<p) in the lattice models described by the bosonic theory. 



A. XX model 



As explained in the appendix, the XX model is mapped into the lattice fermion problem ([B]). In the continuum 
limit, this becomes a massless Dirac fermion whose low-energy regime is exactly formulated, via the bosonization 
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I 1 1 1 1 1 

0.1 0.2 0.3 0.4 0.5 

x = l/N 

FIG. 6: for a right-left particle-hole excitation in the XX model with N = 120 sites. Two values of the external field J are shown: 

for J = and J = 1/2 the filling fractions are 1/2 and 1/3 respectively. Continuous lines (indistinguishable from the numerical data for 
n = 2) are the CFT prediction. 



procedure, as a free bosonic theory (pl5")). According to this map, e 1 ^ (e l< ^) corresponds in the fermionic language to 
the creation of a fermion (hole) (c.f. table [I]). 

Let us illustrate the law Fy _ (x) — 1 for three different states of the finite-size XX problem. Consider |V^,a) with 
(a, a) = (1,-1). It is the Umklapp excitation (c.f. table QJ which, as said in section UH presents a vanishing excess 
of entropy. The same applies for the state (a, a) — (-1,0), which corresponds to the (a) excitation; to (a, a) = (1,1) 
which corresponds to the (b) state (c.f. figure [T]), and, in general, to all the low-energy compact excitations. 

We now focus in the state \idcf>) with total momentum P = 2ir/N. In the XX/free fermion model it corresponds 
to the particle-hole excitation (c.f. table UJ . In figure [5] we present the quantity F-%^ for lattices with different sites 
and filling fractions, together with the CFT predictions (|57|) . (|58p. A further prediction is shown in figure [5] for the 
excitation dcj)d(f>, which corresponds in the lattice to a particle-hole excitation on the top of both the right and left 
Fermi points (c.f. table [I]). The quantity F^^l is compared with the CFT prediction, which is the square of (|57I58|) . 
We observe a very good agreement. 

Finally, in reference [l[ it was shown how, in the x « 1, n — > 1 limits, (|4]) results in the law (|T0|) . which was 
tested in the XX model in figure [21 Interestingly enough, and as can be seen in the figure, the law (I10|) works also 
for non-primary states. 

B. XX Z and excluded-XX models 

We define the XX Z spin chain: 



H XXZ =--Y J [ofo* +1 + °y j+1 + A^f +1 ] , (59) 

3=1 

where PBCs are assumed (for A = one gets the XX model). This model is integrable [23[ and gapless for 
— 1 < A < 1. In the continuum limit it is described by the aforementioned bosonic CFT compactified in a circle 
of radius K = 7r/(2cos _1 (A)). We had access to the exact entanglement entropy of this model through its exact 
diagonalization. 

We first consider the vertex states \V a ) in the XXZ model. The studied excitations were the following: the (a) and 
(fr)-types defined in table HI the Umklapp state and the states \V a ) with a = 1,2,3 (c.f. lj). We had access to the 
exact wave-function of the desired states (and hence to the entropy) through exact diagonalization. We verified that, 
as predicted by CFT, the n-Reny entropy coincides (up to oscillations) with the ground state entropy for n = 2,3, in 
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systems with several critical values of A, filling fractions and sizes up to N = 30 sites. 

Secondly, we studied the state \id<j>) which is the lowest energy one in the sector with total momentum P = 2-k/N . 
The function F' 2,3 ) is compared with the CFT prediction in figure [7]- (a). 

The excluded-XXZ model [24], [25[ is an exactly integrable extension of the standard XXZ chain. In this model 
the up spins (cr z -basis) have an effective size t + 1 (t = 0, 1, 2, . . .) in units of lattice spacing, which means that up 
spins are not allowed at distances smaller than t + 1 sites. The model is defined by the Hamiltonian 



H = ^ E P * ( a j a j+i + + ^)°U)Pt , (60) 

i=i 

where A is the anisotropy and P± projects out states in which pairs of up spins are at distance smaller than t + 1 
sites. For t = one recovers the standard XX Z model (JSHJ) where we only exclude up spins at the same site. Like 
the XXZ model there is a U(l) symmetry, due to the conservation of the number of up spins. The model is critical 
for —1 < A < 1, and is ruled by a Luttinger liquid CFT (c = 1) whose Luttinger parameter K depends on the value 
of t, A and the density of up spins. For A = 0, i. e., the excluded XX chain [24[ K is known analytically, namely 
K = (1 — tp) 2 where p is density of up spins (0 < p < j^)- 

For t = 1,2 and A = we find similar results as for the XXZ model. Figure 0(b) shows some numerical results 
for the vertex and the particle-hole excitation. 




FIG. 7: Entropy ratio i^gi, r = 2,3, for the particle-hole state of the XXZ model for different sizes and values of the anisotropy A in 
the critical region (left) and for the excluded-XX model for t = 1,2 and several filling fractions (right). On inset the Umklapp excitation 
is studied instead, which confirms the law Q51I . 



C. Universality of Renyi entropy oscillations 

Specially important are the models described by Luttinger Liquid field theories, which are CFTs with central 
charge equal to one. Their ground states obey hence the law ([2]) for the Renyi entropy with c = 1. On the other 
hand, finite-size ground states of these models were shown in reference [9j to obey ([2]) up to some oscillations whose 
amplitude turns out to be related to the Luttinger Liquid parameter which depends on the microscopic couplings of 
the Hamiltonian. In this subsection we investigate the oscillations of some elementary excitations of the XX and 
XXZ models. 

For the ground state, the deviation with respect to the CFT prediction ((2]), S^ FT , takes the universal form [9I [26l [27i ] : 



S n (l) - S% FT (£) - f n T n (£/N) cos(2k F £) \D(£,N) sin k F \^ n ' K \ (61) 

where kp the Fermi momentum, D(£, N) = ^ sin7r^/7V is the chord distance, f n is a non-universal constant and T n 
is a universal function depending in general on the parity of N. The exponent p(n, K) = —2K/n is a function of the 
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Renyi index and of the Luttinger Liquid parameter, K. In the models studied in [1, [H, H3; 110 oscillations are found 
for the von Neumann entropy (i.e. for n = 1 in (|5T1) ). We also observe this fact for all the excitations considered 
in this article. In the present section we will provide numerical evidence that the same behaviour (|61j) holds for 
low-energy excitations, the function T n depends in a universal way on the given excitation and the exponent p(n, K) 
is given by —2K/n, at least for some of the excitations considered. 




FIG. 8: (a) Function /„ T n corresponding to the particle-hole excitation of the XX model with four different values of the field J. n = 2, 3 
cases are shown (the former being the upper curve at the half of the chain), (b) idem for the right-left particle- hole excitation. Colors are 
as in figure [6] 



For the XX in an external z-field of strength J (c.f. equation (|A1|) . 7 = 0), the law (joT| holds for all the 
excitations considered. For this model K = 1 regardless the value of J. The iV-dependence of S n — S^ FT at 
half of the chain shows a power law behaviour compatible with p(n, 1) = — 2/n for n = 2,3. Further numerical 
evidence can be found in figure [5] for the non-chiral and chiral particle-hole excitations: the same function f n !Fn 
can be found for common excitations in systems with different values of J and filling fractions, showing the valid- 
ity of equation (pTT) for K = 1, the function T n depending only on the excitation, and withp(n, 1) = —2/n for n = 2, 3. 




x = e/N 

(a) 



x = e/N 
(b) 



FIG. 9: Functions fnFn describing the amplitude of the Renyi oscillations for several excitations in the XX and XXZ (A = —1/2) 
models with N = 240, 500 and 28 respectively, (a) and (b) correspond to the (a) and (b) excitations defined in table U The universality 
of T n is apparent, being this function equal with the ground state case, which is shown for comparison. 



On the other hand, for vertex excitations one observes the validity of the scaling (|6Tj) with p(n, K) also for K 7^ 1. 
A first illustration of this fact can be found in the Umklapp excitation: as mentioned before, the state |Vi,-i) of 
the XX model has exactly the same Renyi entropy as the ground state in the lattice. Moreover, the same state in 
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the XX Z model exhibits the same entropy of the ground state up to corrections of the order of 10~ 3 (i.e., the same 
oscillations) for systems with 30 sites [l| (see also figure [7J- (a)). This implies that, for both XX and XX Z models, 
J- n is common to the ground state and the |Vi,_i) state. Further numerical evidence is shown in figure |9]- a, b, where 
the (a) and (b) states in the XX case are compared with the XXZ case with A = —1/2. The function F n for vertex 
excitations turns out to be the same as for the ground state (apparently a constant function, c.f. figure|9l in agreement 
with [9(). We obtain analog results for other values of A and for n = 3, with p(n, K) — —2K/n as in (l6Tj) . The case 
of the particle-hole excitation for A ^ is more controversial up to the sizes we have investigated, and will be studied 
in a future work which is in preparation [28j |. 

V. FERMIONIC THEORY AND ENTANGLEMENT IN THE CRITICAL ISING MODEL 

The second CFT we examine is the critical Ising model, whose action is given by 

A bP] = ^ / dz d * [W> + (62) 

where ip is a free Majorana fermion. The solution to the equations of motion being ip(z, z) = ip(z) + tp(z). This theory 
is conformal, with c = 1/2 and with primary operators er and ip having the conformal weights: 

(h,h) <r = (±±), (M)* = (^,0). (63) 



Let us first study the function Fj 1 '. One can compute [22 1 the correlators of the field a by bosonization. Two 
Majorana fermions (labelled 1 and 2) can be combined into a Dirac Fermion ^(z) = 2~ 1 / 2 (-0i(z) + iifo^z)). \t is in 
this way a (1/2,0) primary field with two-point correlator: 

(^(z 1 )^(z 2 )) = M*iM*2)) = — ■ (64) 

The CFT describing the Dirac fermion has as central charge twice the Majorana fermion central charge, i.e., c = 1. 
The Dirac theory can be mapped into a free bosonic theory with action (|^5]l in such a way that the spectrum of both 
theories coincide and that to each operator in the Dirac theory there corresponds an operator in the bosonic theory 
with common algebra and correlators. For the field "J, this mapping is: 

$(z)=e'*W, *(z) = e^ (2) 
and for a (as can be checked from the matching of the conformal weights): 



0"! 



(z, z)a 2 (z, z) = 2 1/2 cos(tp(z, z)/2). 



As described in [22|, the following trick is used for the computation of a chain of a correlators. Two copies of a 
operator chains are constructed 

_. . _.. , . . / (p(z,z) ip(w,w) 
(<7i [Z, Z)<J2(Z, Z)(Ti (W, W)<J2\W, w)) = {<J(z, z)<j(w 7 w)) — 2 ( cos — - — cos 



Using the vertex operator correlator (j48|) one arrives to: 

(a(z, z)a(w, w)) 2 = \z — w| -1 ^ 2 , 



and, for a set of 2n fields, to: 



(^(z 1 Z 1 ) . . .a(z2 n Z2n)) 2 = 2 " ^ Yl " 

{a J= ±l} i<3 



iaj /2 
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where the sum is subject to: Ylj a j = 0- I n the cylinder w = — ilnz this is: 

(cr(wiwi) . . . a-(w 2 nW2n)) 2 = 2~ n | sin ? 

{q j= ±1} i<j 

With this information one can construct the function F„ , which reads: 



(65) 



= |shu>x)|"/ 2 

{a 3 =±l} 

n — 1 n — q 

nn 



7r(q — x) 
sin 



n— 1 n— <j 

nn 

9=0 4=1 

'2j°2(j + i;)-l 



. 7r(o + IE) 

sin 



*2j-l°2(j + <,) 



sin — 

n 



»2j°2(j + g) 



1. 



(66) 



9=1 4=1 

This function can be proved to be constantly equal to one for all values of n, in the same way that this fact was 
proved for Fy . The square of equation (|66| can be written as the product of two terms: the function |sin(7ra;)|" 
on the one hand, which presents a zero of order n for each igZ, and the sum on the permutations squared, which 
presents poles in x G Z of order n. Indeed, for each integer whose rest of the division by n is p > 0, there is a pole of 
order n — p in the q = p product in the second line of (|66p . and a pole of order p in the q = n — p product in the first 
line, and this happens since, for a given p, the permutation {aj = (— 1) J } satisfies 2j=i a 2j a 2(j+p)-i = n — p and 
Sj=i a 2j-io.2{j+n-p) — P- We hence find the CFT prediction that the entropy of a excitations coincides again with 
the ground state entropy. 

We repeat the game for the thermal field s(z,z) = iip(z)ip(z), with conformal weights (h,h) = (|, |). For the 
computation of F^ we use the 2n-point correlator of e's (cf. (22[, p. 444), which can be inferred from the ip correlator 
(El: 



(e(zi,Zi) . . .e(z2n,Z2n)} 
{lp{zi) . . . 1p(z 2n ))(lp(zi) ■ ■ ■ 1p(z 2 n)} 



Pf 


1 






z i3 _ 


l<i,4<2n 



where Pf means the Pfaffian: 



(67) 



1 n 

pf(A)= 2^ n A ^i-iu 



(2i)- 



In the cylinder w — — zlnz and taking the 2n coordinates ro,j as in one arrives to: 



(68) 



2n 



4=1 





1 




^Pf 






l<i,4<2n/ 



(69) 



and using that the Pfaffian of an antisymmetric matrix is the square root of its determinant, one finds that Ft\x) 
coincides with the value computed in section |IV] for the field id(j> (cf. equation ([56]) ): 



F^(x) = F%l(x). 



(70) 



Finally, using (|67|). one obtains 



Fl n) (x) = ^F%l(x). 



(71) 
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A. Ising model on the lattice 

Consider the critical Ising model on a transverse field (ITF) 



1 



JV 



3=1 



J a* 



(72) 



The model is integrable and critical for J = 1. The critical model is described in the continuum limit by the free 
Majorana fermion (p)2"]) (22[, as can be seen from the spectra of its free-fermionic formulation (jAllI) . We identified 
the states \a), |e), \ip) in the finite size Ising model through their finite-size scaling of energy and momentum and for 
these states we computed the entanglement entropies and the quantity via the method described in [2!| and in 
the appendix. 



i?(2) 



PI □ 



N = 242 
JV = 202 
JV = 162 
JV = 122 
JV = 82 
JV = 42 



H n 



° t V% B ; S ■ ■ ■ . ■ BBBDDbdbbbi; 

° <> o . ^OOGO<J 



0.1 



x = l/N 

FIG. 10: The function seems to converge to the CFT prediction Fcr(x) = 1 in the large JV-limit of the critical Ising model. The 

inset shows the JV -1 ' 2 -dependence of F^ 2 '(l/2). 



The excitation \o) corresponds in the lattice to the lowest energy state in the parity -1 sector. Figure fTOl shows the 

(2) 

quantity F# for systems of several sizes up to N = 242. One observes that, although the iV-scaling is slower than in 
the bosonic case, the curve seems to converge to the CFT prediction for large N (see N scaling of Fe(l/2) in 
the inset). For Fj 3) we obtain analog results. 



On the other hand, a comparison between the numerical values for the quantities F 



(n) 



2, 3 and the laws 



(I70I71P is shown in figure [TT] Once more, the agreement between CFT and finite-site lattices is remarkable. 



VI. DESCENDANT FIELDS 



The results presented so far give a strong numerical support to the validity of equation (|43l) for the entropy ratio 
F^ when T is a primary field. A natural question is whether that formula also applies when T is not a primary 
field but a descendant. We have not yet found a definite answer to this question but there are some observations 
that may eventually lead to a solution (for details see [2^]). Let us first recall that in the field theory construction 
of the reduced density matrix, px, one inserts into the path integral an operator at r = — oo, and another one at 
t = oo, in order to create the incoming and outgoing states associated to T and Y 1 * . These operators are inserted 
in a cylindrical geometry and hence the corresponding operators must be the conformal transformed of those defined 
on the complex plane. For primary fields, the latter conformal transformation is rather simple but for descendant 
fields it generates additional fields in the conformal tower. A typical example is the energy-momentum tensor, T, 
whose conformal transformation, from the complex plane to the cylinder, involves a constant term proportional to 
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FIG. 11: n = 2 and 3-Renyi entropy ratio for the states \ip), \e) of a N = 200 critical Ising, compared with the CFT predictions for these 
fields (continuous curves). 



the Schwarzian derivative. Hence in this case the operators to be inserted at r = ±oo involve T plus a constant term. 
However this constant term cancels out in the conformal transformation (I29|) from the cylinder to the uniformizing 

plane z, so that the final expression for F^ n) is given by equation (|39[) with T replaced by T in the 2n-point correlator 
on the complex plane. The same result can be obtained for the descendant field d 2 <j) of the c = 1 theory, and we 
expect this result to hold in general. It would thus seem that equation (f3"9")) must also be valid for descendant fields. 
However the numerical results give only a partial confirmation of this conjecture, which we now explain in more detail. 

Let us first consider the entropy of descendant fields in the sector with (h,h) = (2,0) of the bosonic CFT. This 
space is the first (two-dimensional) degenerated sector of the identity tower in the bosonic CFT, and it is expanded 
by the fields : (dip) 2 : and : d 2 ip : (the former being proportional to the stress-energy tensor of the theory, T). In the 
XX and XXZ models this sector has total energy and momentum E = P = 4n/N. 

In reference [l[ it was noticed that the entropy of degenerate sectors could depend on the particular state in the 
sector. We will provide a example of this fact. Let |1), |2) denote a basis of the (h, h) — (2, 0) sector in a finite-size 
realization of the XX chain. We computed the entropy S2 for a general eigenvector of the XX model in the mentioned 
sector: 

\v) = cosa|l) +e l/3 sina|2) (73) 

via exact diagonalization. Figure [U shows FrJ for several values of a and f3, together with the entropy of the state 
(1 : 2) computed with the correlation matrix method exposed in the appendix (circles). The states (1 : 2) and (2 : 1) 
constitute a basis of the sector, an both have the same entropy. Different grey levels in fig|T2] denote different values 
of a. There is a certain linear combination \v) for which the entropy coincides with the entropy of the state (1 : 2), 
which turns out to have the lowest value. The CFT prediction for T, obtained using equation (f3"9"| . is also reported 
in figure [T2l and it seems to coincide with one of the states \v). 

As another example, let us consider the descendant state (1 : 2), whose field theory counterpart must be a linear 
combination Y of the operators 

Y = cos a d 2 <j) + e ip sin a T (74) 

The values of a and /3 can be chosen to fit the numerical values of t° the value of Fy obtained using equation 
(|39[) . The result shown in fig |13l shows a remarkable agreement, which however is not as good as in the case of primary 
fields. Another example is the entropy of the energy-momentum tensor of the Ising model, where the CFT prediction 
given by equation (|39[) disagrees strongly with the numerical results. All this shows that the general application of 
this formula for descendants fields deserves further clarification. 
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FIG. 12: n = 2 entropy ratio of several states in the (fc, fc) = (2, 0) sector of the XX model with N = 16. 



0.9 



0.8 

i?(2) 

0.7 



0.6 



0.5 

0.1 0.2 0.3 0.4 0.5 

FIG. 13: n = 2 entropy ratio corresponding to the (1 : 2) excitation of the XX model with N = 400 (points). The continuous curve is 
the a fit with the function Fy with Y given by 174> . 



VII. CONCLUSIONS AND PERSPECTIVES 

The main findings reported in this article are summarized in what follows. The universality of the ground state 
Renyi entropy of entanglement at criticality is generalized to low-energy excitations represented by primary fields in 
CFT, as stated in reference [l|, the n-th Renyi entropy being related to the 2n-point correlator of the corresponding 
field. In this work, this result is derived in the more general framework of the Calabrese and Cardy article of 2004. The 
result is found to reproduce correctly the entropy of primary states of several models: the XX, XX Z and excluded- 
XX Z models in the c = 1 universality class, and the Ising model in the c = 1/2 universality class. Predictions for the 
low-£ behaviour of the von Neumann entropy are also numerically verified for bosons. We present numerical evidence 
of the universality of the parity effect of the Renyi entropy also for excited states, although further investigation is 
needed. 

Finally, we have studied the Renyi entropies of some descendant states in the XX and the Ising model, finding a 
partial sucess of the analytic approach which indicates that further investigation is required to fully understand the 





20 



entropies of all the low energy states in critical lattice systems. 
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Appendix A: Entanglement in the XY model via correlation matrices 

We consider the XY model for N spin-1/2 particles with periodic boundary conditions (PBC) in an external z 
field. 

-1 N 1 

HxY = ~2 E o [(1 + 7) ^ fT J+ 1 + (1 - + Ja l ( A1 ) 

where a a are the Pauli matrices: 

-(?;)■ -{iiy 

This generalizes the Ising (ITF) (7 = 1) and XX (7 = 0, J = 0) models. Defining the raising and lowering spin 
operators a ± = \{a x ± ia v ) we have: 

-1 N 

HxY = — E + 7«+i + h.c] + Ja| 

i=i 

A Jordan- Wigner transformation is now performed to express (| Al[) in terms of true fermions {c^, c„} = <5 mj „. They 
are defined as: 



n ^ ] °™ ( A2 ) 

and they satisfy 

cjyci = -(-l)"V+ar (A3) 
where = N — uf is the number of down spins and 

N z 1 1 
SJ + 1 



is the number of up spins (s| being the eignevalue of and also the number of fermions. Equation (|A3|) determines 
the boundary conditions in the fermionic formulation of (|Alj) and the Hamiltonian reads (cn+i = (— l)" J 'Ci): 



1 W 

" 2 E [ c i c J+i +7c)ct +1 + h.c. - 2Jct Cj - 
3=1 



A Fourier transform is now performed. The N fermionic modes, {dt,df} = <5j-j/ are defined: 
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1 N 

^■ = 7iE e ^ rac - ( A5 ) 



N 

m—1 



(where kj = 2tt j/N is the momentum associated with index j), in such a way that: 



H XY = ^2 A i d \ A i + iB i d ] d -j + d o d i ■ ; Mii 



where we defined Aj, Bj 



7 

Aj = cos fcj — J, = — sin kj , (A7) 

and where the set of momentum indexes f2 is such that the resulting allowed momenta kj, j £ CI, are those specified 
later in equation IA131 The Hamiltonian (|A6I) describes interacting fermions and can be diagonalized throught a 
Bogolubov transformation. Defining the fermionic modes (v3 +i>| = 1): 



b] = Uj d] + ivj d-j 
b-j = Uj d-j + ivj dj, (A8) 

and imposing that the Hamiltonian for the j, —j modes has the diagonal form (Ej bjbj — Ej bjbj)/2, i.e., that: 



j j 
we obtain the Bogolubov equations to be satisfied by Uj , Vj : 



and (|A9t with 



-UjVj = Bj/Ej (A10) 



E] = (cos kj - J) + 7 2 sin 2 kj. (All) 

There is an ambiguity in the sign of Ej in (jAlip . Each j, — j being diagonalized independently, one can choose the 
sign of Ej independently, Ej — E_j — \Ej\sj. With the election Sj = —1, the ground state of Hxx corresponds to 
the ground state of the free fermion system (|A9I) at half filling, as said in section [TTI (even N is supposed). 

In conclusion, via the Bogolubov transformation (|A8IA9IA10[) . we have expressed the XY problem in terms of a 
free fermion problem, in such a way that the eigenstates of Hxy in this free fermionic formulation are: 

n «4i°> ( ai2 ) 

and that the spectra of Hxy is X^eK Ej ~ K, being |0) the fermion vaccum and K C CI is one of the a set of np 
integers or half-integers characterizing the state, such that kj = 2tt j/N is the momentum corresponding to index 
j £ K. The set of allowed momenta depends on the boundary conditions of the fermionic problem, which are in their 
turn determined by the parity of (IA3|) : 



N-n F even 
(APBCs) 



N even: N(j) = ±tt, ±3tt, . . . , ±(N - 1)tt 
N odd: Ncj) = ±tt, ±3tt, . . . , ±(N - 2)tt, Ntt 



VbCs^ { N even: N<f) = °' ±2 ^' ±4 ^' ■ • ■ ' ±{N " 2)?r ' Nn ( A13) 

1 N odd: N(f> — 0, ±27r, ±47r, . . . , ±(N — 1)tt 
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1. Correlation matrix in terms of Majorana fermions 

The 2N spatial Majorana modes c r are denned: 

C2m-1 + 1C2m 



C2m-1 
C2m 



(A14) 



they satisfy Majorana anticommutation rules: {c r ,c s } — 25 rs . We want to construct the correlation matrix (c r c s )n 
for a given state defined by K. We first express the c's in terms of the d's, then the c's in terms of 6's using the inverse 
Bogolubov transformation (|A8[) . Finally, the set K determines the occupancies of the free fermions: (bjbji)n — Sj j> 
if h j' G K, and zero otherwise, and (bjbj, ) = (bjbf) = 0, and one gets for (c r c s ): 



with r e M2 



(A15) 



/ n Hi 
n_i n 



Vn. 



N+l 



n 



JV-l 



n / 



(A16) 



and where: 



/ o (1) o (2) 

TT _ I yro ym 

Um ~ „(3) „(1) 



(A17) 



7 (1) = ± 



Ask 



e + > e 

fe^K 



7 (2) _ 



E 



-i<t> k m 



7 (3) 



iV 



iV 



fe^K 



.feeK 



e~^ m (u k + iv k f - Yl ^ fcm K - ^ fe ) 2 



(A18) 



LfceK fc^K 

From the antisymmetry of (c r c s ) one can see that gffl — —g^m and g 1 ^ = —g^l} n - For the ground state, it is 
K = — K and the same for the complementary of K (this is not true for a general state). One now uses this condition: 
J2j S-j m the expression for g m , and also the fact that, for the ground state, Ej = — \Ej\ if j E K, E k = \E k \ 
otherwise, and we obtain: 



1 E ^ 



2iB 



N 



all 



(A19) 



which in the continuum limit becomes equation (3.47) in [30(, of which (|A15|) is a generalization for arbitrary states 
characterized by K and for finite-size systems. 
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2. Computation of entanglement 

We are interested in the entanglement between a spatial partition containing the first £ sites of the system, and the 
rest of the system. Being ci-m, , C2 m -i related to the m-th site, the reduced density matrix of the subsystem defined 
by 1 < £ < N is encoded in a 2^-block of matrix T, Y e € M 2 e(C): 



( n Hi 
n_i n 

Vn_ £+1 ... 



n, 



(A20) 



n / 



In particular, the entanglement is computed in the following way. As ((c r c s )) rs , Yg is an antisymmetric matrix whose 



eigenvalues are complex coniugated pure imaginary numbers ±iv r 
matrix can be transformed into a block-diagonal form: 



being real and depending on I. Such a diagonal 




(A21) 

Now suppose that a m is the basis in which this happens, i. e., (a r a s ) = S rs + i(Yi) rs . It is easy to see that the 
corresponding true fermions a m = a,2 m -i + ia>2m are in a product of uncorrelated states: 



(a m a n ) = 0, (a m a n ) = 5 mn — (A22) 

This product being uncorrelated, the reduced correlation matrix of the ^-block is a product ® m Q m of correlation 
matrices of single-site blocks Qj = pj a]|0)(0|aj + (1 — pj)\0)(0\, pj = (1 + Vj)j2. The entanglement entropy S(£) of 
the ^-lengthed block is, hence, the sum of the entropies of the £> m 's: 



<?W = £if 2 (^p). (A23) 

being H2(x) = —x In x — (1 — x) ln(l — x). The strategy is hence to construct the matrix (|A20[) numerically for a given 
state caracterized by the set K. By diagonalyzing Yi we obtain v r and finally the entropy via (IA23p . In the same 
way, the Renyi entropy can be computed: 



l 

5„W=^ln[p/ l + (l- Pj )"]. (A24) 

3=1 

The total time cost of the operation being 0(£ 3 ). 
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